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t^J- \ the Nambu-Jona-Lasinio model in its bosonized form. We derive current quark mass 



We explore the influence of the current quark mass on observables in the low 
energy regime of hadronic interactions within a renormalization group analysis of 



expansions for the pion decay constant and the pion mass, and we recover analytically 

>: 

• the universal logarithmic dependence. A numerical solution of the renormalization 

group flow equations enables us to determine effective low energy constants from the 



model. We find values consistent with the phenomenological estimates used in chiral 
perturbation theory. 



P-i 

qj\ I. INTRODUCTION 



Chiral Perturbation Theory (%PT) is an effective field theory which describes the low energy 



limit of the physics of the strong interaction where its theoretical description, Quantum 
Chromodynamics (QCD), is entirely non-perturbative. It is based on the expansion of an 
effective Lagrangian in terms of small external momenta and the pion mass. A finite pion 
mass arises from the small explicit breaking of chiral symmetry. As an effective field theory, 
xPT is defined for scales much smaller than the chiral symmetry breaking scale Av ~ 1 GeV, 
and it is renormalizable order by order in the small momentum expansion Q, Q, [J- The 
physics at low momenta is completely dominated by the light pions. Other fluctuations are 
separated from the light Goldstone bosons by a mass gap and therefore suppressed. Due to 
the constraints of chiral symmetry, all quantum corrections are due to pions and depend on 
the pion mass parameter. However, going beyond leading order in the chiral expansion, the 
number of undetermined coupling constants in the chiral Lagrangian increases rapidly. 
The chiral expansion in small pion masses and momenta gives rise to terms that are non- 
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analytic in the pion mass, i.e logarithmic corrections appear at the one-loop level. In leading 
order, these are universal corrections and entirely determined by the symmetry. Therefore, 
they should be reproduced by any theory that describes the low energy domain within a 
suitable approximation. On the lattice, the appearance of these chiral logarithms has not yet 
been unambiguously demonstrated, but recent results remain compatible with the existence 
of such corrections 4J . It is our purpose to show that the logarithmic terms are recovered in 
the pion sector of the bosonized Nambu-Jona-Lasinio (NJL) model 0,0], and to determine 
them in a Renormalization Group (RG) approach. 

An effective field theory can be considered as the limit of a more fundamental theory at 
low energy. Only fields relevant in this energy range are retained, while all other degrees 
of freedom are integrated out. The effects of these higher modes enter into low energy 
couplings (LECs) which multiply the remaining operators in the low energy theory. These 



LEC 8 have to b e obtained either fr o,„ exponent tj, o r ov non-pe—ve .nethods fr0m 
the underlying gauge dynamics of QCD. Because of the strong gauge coupling at low energies 
and the emergence of collective hadronic degrees of freedom this remains a difficult task. 
Recent progress in a non-perturbative description has been made using Dyson-Schwinger 
Equations (DSE) [tUsH^, or RG flow equations 10, 11, 12, [l4|. Numerical lattice simu- 
lations are limited by the computational effort and require extrapolations to large volumes, 
small quark masses, and small lattice spacing. The extrapolation to small quark mass is 
actively pursued |4| and an important topic for yPT. Finite volume effects have been ad- 



dressed in xPT as well, both for the nucleon |15j| and for the pion 



the framework of the quark meson model with RG methods 
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and in 



Due to the importance of the low energy couplings on one hand, and the difficulty of ob- 



taining them from QCD on the other, there 



models 
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ras been a strong interest in phenomenological 
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351 ] . RG methods are uniquely suited 



to interpolate between the low energy effective theory and QCD [33]. 

There are bottom-up or top-down models. In the bottom-up approach, the starting point 



is X 
ref. 



'T. Additional vector, axial-vector and scalar meson resonances may be included. In 



24j , it was argued that the coupling constants of the chiral Lagrangian in next-to-leadin 



order could be accounted for mainly by contributions of vector meson resonances. In 
it was found that in a linear meson model the low energy couplings are dominated by the 
exchange of scalar mesons. Such descriptions rely entirely on a hadronic picture. On the 
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other hand, the top-down approach starts from models with quark degrees of freedom, which 
generate mesonic degrees of freedom by quark dynamics. These approaches are essentially 
based on the NJL model , which includes a fermionic self-interaction. Meson fields arise 

rom a bosonization of the quark interaction with a Hubbard-Stratonovich transformation 

~l 

3(t| . Pions appear explicitly as the Goldstone bosons of spontaneous breaking of the chiral 



symmetry through a quark condensate. For a review of this class of models, see e.g. 
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321 ] . More recently, 



A comparison with xPT has been undertaken e.g. in 0,0,0, 
a third approach based on a weak form of the ADS/CFT correspondence J37J shows a 
remarkable agreement with hadr qn p henomenology 3sl |. 

In the RG treatment 0,0, 41, Q] of NJL-based models, chiral symmetry is dynamically 
broken, and the pions emerge naturally in the RG flow of the bosonized model. In order to 
capture the low energy dynamics dominated by the pions, the meson dynamics has to be 
included in the RG flow. This goes beyond the leading order of the often employed large- N c 



approximation 



n 



Only the inclusion of such 1/N C corrections [42] makes it possible 



to derive the leading-order logarithmic corrections in the quark mass expansion of low energy 
observables. In particular, our results show that the wave function renormalization for 
the derivative term of the mesonic fields is essential, since it enters into the expression for 
the pion decay constant /„-, the coupling constant for the lowest-order derivative term in 
xPT (compare j^). 

In a comparison of RG results to %PT, one would like to calculate the same effective cou- 
plings. Already the lowest order xPT Lagrangian contains information about arbitrary pion 
n-point functions, which makes an order-by-order comparison with the results of RG meth- 
ods difficult, since the effective RG action is expanded in n-point functions. Alternatively, 
one can proceed in the same way as in the comparison with experimental results. By calcu- 
lating observables and identifying the influence of the low energy constants one can compare 
effective values for these couplings. In this way one does not actually have to project the 
RG flow onto the operators used in the chiral expansion. We will employ this method in 
the present paper. Our comparison between xPT and the NJL model will be limited to the 
two flavor case and of course to only those observables that appear in both theories. In the 
specific case of the bosonized NJL model, which we define at a large UV scale, there are 
several free parameters. The UV parameters are fixed by requiring that the model should 
reproduce the physical values for the pion decay constant f n and the pion mass m n , after 
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all quantum fluctuations have been integrated out. While the actual values of f n and m 7T 
are input parameters to the model and not predicted, their predicted dependence on the 
current quark mass tests whether the model as a low energy description of meson physics is 
compatible with other approaches. Since the low energy constants depend on the underlying 
short distance theory, they differentiate whether the model contains the relevant dynamics 
for the emergence of the low energy physics. Such a consistency check can put the model 
on a more sound footing, even without a direct connection to the gauge dynamics of QCD. 
The main part of this paper is organized as follows: After a short review of some results of 
chiral perturbation theory at one-loop order in section |Hj we will present the RG approach 
to the bosonized NJL model and introduce the RG flow equations in section IIHI Analytic 
results for the non-analytic dependence of the pion mass and the pion decay constant on 
the symmetry breaking parameter are given in section IIVI Numerical results on the same 
quantities as function of the current quark mass are given in section We will demonstrate 
that the results from the NJL model for the couplings that appear in the pion mass and the 
pion decay constant are compatible with the values used in calculations. We present 
a summary and our conclusions in section IVI1 Details of the derivation of the RG flow 
equations can be found in the appendix. 



II. RESULTS FROM CHIRAL PERTURBATION THEORY 

In this section we give a short summary of some principles and results of at the one- 
loop level following the work of Gasser and Leutwyler and the recent more pedagogic 
article Starting point is the Lagrangian £ x pt of the non-linear sigma model in the 

0(4) chiral field U = (U°,U l ) of unit length expanded up to fourth order in momentum. 
We give this Lagrangian with the 10 phenomenological coupling constants (k, hi) here for 
later comparison with the RG-results 

+ h {V^U T V V U){V,U T V V U) + h (x T u) 2 + h (v^ T v M f>) 
+ h (U T F^F,„U) + k (VUTF^VU) + l 7 (x T U? 

+ h X T X + h 2 trF^ v + h 3 x T x, (1) 
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with 

V M f> = d^U° + aft* 

V M C7* = fyt/* + e ikl v k Jj l - a^U 

X = 2B(s°,p*), 

X = 2B(p°,-s i ) 
F^U = (V^V, - V„V„) £>, 

where a* (x), s(x) and are external vector, axialvector, scalar and pseudoscalar 

fields. Without isospin breaking the constants I7 and are zero. 

The constant f n o is the pion decay constant in the chiral limit, which in the literature of 
chiral perturbation theory usually is called F, and B parametrizes the chiral condensate in 
the chiral limit, 

1 (gg) = 1 (uu + d d) Q 

In the two flavor case treated in this work, we consider the masses of the two lightest flavors 
to be equal. In order to take the quark mass term into account one expands the scalar 
external field around s = m c , where m c is the average current quark mass: 

m c = ^{m u + m d ) (3) 

One obtains the Gell-Mann-Oakes-Renner relation, where the parameter M is the pion 
mass in lowest order in the current quark mass 



M 2 = 2m c B = -m c M^. (4) 

/ttO 



The non-linear sigma model is not renormalizable, but it is possible to make the theory 
finite at every loop order by introducing appropriate counter terms. In chiral perturbation 
theory this renormalization is usually done by dimensional regularization. Up to one-loop 
order, one needs the above ten low energy constants and hi which depend logarithmically 
on the renormalization scale /1, 



K = ^{h l + \og^) 1 = 1,2,3. (6) 
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The low energy constants U are independent of the renormalization scale and are related to 
physical observables. The constants hi have no direct physical relevance, but follow from 
the renormalization procedure. To go to higher loop order, an increasingly larger number of 
constants is needed. By computing Greens functions from chiral perturbation theory in the 
limit of very small external momenta, one derives physical observables which are expanded 
in M 2 . As shown above, M 2 is related to the current quark mass m c . For future reference, 
we quote the chiral expansions for the pion decay constant, the chiral quark condensate and 
the pion mass, 



We will compare our RG results with these expansions. 

The expansions are characterized by the so-called chiral logarithms and by higher order 

terms in the expansion parameter M 2 . It is possible to combine the chiral logarithms with 

the linear term in M 2 (cf. equations (jSJ) and (jHJ)), so that the constants which appear do 

not depend on the renormalization scale [i any more, but only on M 2 , 

( M 2 \ 

U = /*o 1 + TTrnrk + 0(M 4 ) (10) 



16vr 2 / 2 



o 



/ /if 2 \ 

(w> = <??>«, (i + ^jr^ - r ») + ) ( n ) 

< = M 2 (l- -^h + 0(M*)) . (12) 



V 32tt 2 / 2 , 
The renormalization scale-independent constants U are parametrized in terms of the con 
stants A 2 



i ; 1 

,2 



Af T = 0.59+^° GeV , Af T = 1.26 ± 0.14 GeV (14) 

The M 2 -dependence is contained in the logarithms, Ag PT and A;![ PT are scale-independent 

461 ] . while is not determined in chiral perturbation theory since 



constants, as given in 



(qq) is not a physical quantity. 
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III. RENORMALIZATION GROUP FLOW EQUATIONS FOR THE 

BOSONIZED NJL MODEL 

In this section, we give a brief overview of the renormalization group approach to the NJL 
model. The detailed approximation scheme and the derivation of the flow equations are 
presented in l^,!^- Here we only review the main results and give the flow equations that 
are used for the analysis of the chiral low energy expansion in the next section. 
The basic idea of the renormalization group is to describe the dependence of an effective 
action on an infrared cutoff scale, and to follow the evolution of the couplings under a change 
of this cutoff scale. In this way, quantum fluctuations are integrated out in a systematic 
way, and a theory with bare couplings defined at a given UV scale is transformed into an 
effective low energy theory in which all quantum fluctuations with large momenta are already 
included in the couplings 149, 501. 

nn 

Our starting point is the NJL model [5J, |6J with quark fields interacting via a point-like four- 
fermion interaction. We use this model of dynamical chiral symmetry breaking to describe 
the transition from quark fields to hadronic degrees of freedom for scales k ^5 1 GeV. The 
Lagrangian reads in Euclidean spacetime 

£njl = q($ + m c )q - g NJL {(qqf + (qir^q) 2 ) . (15) 

This Lagrangian contains the dominant chirally symmetric four-fermion interaction. Chiral 
symmetry is preserved by the interaction term, but broken explicitly by a finite current quark 
mass m c . By bosonization of the quark fields at the UV scale, we obtain the associated linear 



cr-model 



51 1 , where the chiral symmetry of the quark fields (q, q) is for two quark flavors 



reflected in an 0(4)-symmetry of the meson fields $ = (a, 7?), 

£ uv = 0q + g vv q((T + irif lb )q + t/ uv ($ 2 , a, k vv ) (16) 

with 

rj uv ($ 2 , a, k uv ) = ^(a 2 + vf 2 ) - 5a, (17) 

2 

771 

§ = r ^2L mc , (18) 

The parameter 5 depends on the current quark mass and describes explicit chiral symmetry 
breaking. 
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We sketch the derivation of renormalization group flow equations for the effective action T 
of the bosonized NJL model in the Schwinger proper-time regularization scheme 5^|. The 
effective action is the generating functional of the one-particle irreducible Greens functions 
of quarks and meson fields. The derivative expansion of the effective action up to next-to- 
leading order is given by 



T[q,q,^} = J d A x(^{<b\o) + ^Z£{&){d^ a ){d^ 1 



+G($ 2 ) q(a + irn l5 )q + Z q ^ 2 ) 0qj (19) 

In general, the effective wave-function renormalization factors Z|, b ($ 2 ), Z q (§ 2 ), as well as 
the effective Yukawa coupling G($ 2 ) in this action, are functions of the bosonic fields $ 
and thereby include arbitrary many "elementary" couplings - in full analogy to the effective 
potential U($ 2 , a). In ref. |42( a solution is given for this general case. In this work we only 
treat the effective potential U(& 2 ,a) as field-dependent, while we approximate the other 
couplings by taking only their first, ^-independent term into account, i.e. 

= Z*{al) 6 ab , G($ 2 ) = G(a 2 ), Z q {<$> 2 ) = Z q (a 2 ), (20) 

where <To denotes the vacuum expectation value from the effective potential U. We use 

the same wavefunction renormalization Z<j,(cr 2 ) for all mesonic fields. The flow equations 

for the bosonized NJL model in the proper-time renormalization group scheme are derived 

in ref. j^J. These equations for the individual couplings of the effective action take the 

following general form: 

k dU(&,k) k 6 ( AN c N f 1 I \ 

dk 32tt 2 \k 2 + M 2 k 2 + M 2 k 2 + M 2 J 1 ' 



dZ*(*lk) k 6 ( AN c Nf G 2 2 F 2 

dk ' 16tt 2 $ V (k 2 + M 2 ) 3 Z 2 Z^ (k 2 + Ml) 2 (k 2 + M 2 f 

dG(a 2 ,k) k 6 G 3 ( 6k 2 + 3M 2 + 3M 2 2k 2 + M 2 + M 2 



dk 16tt 2 Z 2 Z$ V (k 2 + M 2 ) 2 (k 2 + Ml) 2 (k 2 + M 2 ) 2 (k 2 + M 2 ) 2 

; dZ q (a 2 ,k) = ( 9 M 2 + Ml _ Qe (M 2 - M 2 ) 2 



(22) 
(23) 



dk 32vr 2 Z q Z^ V (k 2 + M 2 ) 2 (k 2 + M 2 ) 2 (k 2 + M 2 f (k 2 + M 2 ) J 



M*+K „, 2 W -Ml) 2 , (2j) 



+3 _ "«'"' a — — - 2k 



(k 2 + Ml) 2 (k 2 + Ml) 2 (k 2 + Ml) 3 (k 2 + Mlf 

The flow equations describe the evolution of the effective potential U, the Yukawa coupling 
G, and the wave function renormalizations Z$ and Z q under a change of the renormalization 
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f g 2 ($ 2 , 



scale k. The effective masses in the flow equations are given by 

*) = (25) 

(28) 

the effective pion decay constant F 7r ($ 2 ) by 

Fl = (29) 

and the effective four-boson coupling A($ 2 ) by 

2 d 2 U(<5> 2 ) Ml -Ml . , 

Z% (<9$ 2 ) 2 2 2 $ $ 2 ' 1 ' 

They depend on the field <3> 2 and the scale k. Evaluated at the vacuum expectation value 
$ = do, the masses reduce to the physical masses at the scale k 

m q (k) = M q ($ 2 ,k)\^ =aV m^k) = M^ 2 ,k)\^ 2=aV m a (k) = M CT ($ 2 , A;) | 02=ct , . (31) 

A flow equation for the vacuum expectation value (a) = cxq can be derived from the minimum 
condition of the effective potential, 

d 



= 0. (32) 

(7=(TQ ,7T 2 =0 



Details of the derivation can be found in appendix El here we only quote the result 

(33) 



where Uo(<& 2 , k) denotes the symmetric part of the effective potential, without the symmetry 
breaking term proportional to 5, and the prime denotes a derivative in the symmetric variable 
$ 2 . Note that the set of RG flow equations f(|21 )) - (J24jO depends only on the symmetric 
variable $ 2 . Explicit symmetry breaking introduced by the linear term 5a affects only the 
initial condition ao{kuv) in eq. (|HHJ). 

Two generic solutions of the renormalization group flow equations are shown in Fig. Q 
For zero current quark mass m c = 0, spontaneous chiral symmetry breaking sets in at the 
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FIG. 1: The scale-dependent constituent quark mass m q and the meson masses {rria^m^) as func- 
tions of the infrared cutoff parameter k. Shown are results for simplified solutions of the RG flow 
equations with constant wave function normalizations Z$ = Z q = 1 and constant Yukawa coupling 
G. Dashed lines show the chiral limit m c — ► 0, while solid lines give the solutions with explicit 
symmetry breaking. The grey bars indicate the bands in which the masses are approximated as 
constant in the analytical calculation of section IV. 

chiral symmetry breaking scale k x sB and then develops fully in the infrared. This picture is 
characterized by a non-vanishing quark mass m q and a sigma mass m CT 2m q in the infrared. 
Below the chiral symmetry breaking scale, the evolution of the vacuum expectation value a 
quickly freezes. For a finite current quark mass m c , see full lines in Fig. 1. The transition 
from the UV-dynamics to the IR-dynamics is smooth, the meson masses never become zero. 

IV. ANALYTICAL LOW ENERGY EXPANSIONS FROM THE 
RENORMALIZATION GROUP 

In this section, we show that the RG flow equations of the bosonized NJL model generate 
the proper leading order pion contributions to physical observables. Due to the long-range 
fluctuations of the pions, a logarithmic term (chiral log) appears in the quark mass expansion. 
The RG flow equations generate such a term. Their full solution is best suited to sum up 
all long-range fluctuations, as we know from the physics of second order phase transitions. 
Since the RG flow equations go beyond an expansion in powers of the coupling constants, 
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they cannot be easily mapped onto a perturbative expansion. However, we can identify 
terms with the same analytical behavior in the lowest order pion mass, the parameter M, 
and make a direct comparison to xPT results. In particular, the one-loop perturbative result 
can be recovered from RG flow equations by replacing all running masses by constant masses 
in lowest order Treating the small pion mass as constant in the RG flow equations is 
justified only below the chiral symmetry breaking scale k x sB, where the pion exists as a 
low mass and low momentum degree of freedom and pion fluctuations are the dominating 
effects. The chiral symmetry breaking scale k x sB in the RG approach appears as a UV cutoff. 
Although k x sB sets the scale for the logarithmic dependence on the pion mass parameter 
M 2 , additional contributions may modify the value. 

As a consequence of fixing the masses, the flow equations (|21jl and (|22jl decouple from 
equations (|2*3*Jl and (J2*3j) . This can be seen by substituting equations (|2*U|) and (j3*Uj) in equation 



and replacing the ratio ^ by J^) • Therefore the evolution of G and Z q is not relevant 



for the solution of equations (JzlJ) and (j22j) . The system reduces to only two relevant equations 
for the quantities we consider, one for the effective potential U and one for the wave function 
renormalization Z^, 

k d JL = _JL ( ^ N f _ 3 ^ L_] (M) 

dk 32tt 2 \k 2 + m 2 q k 2 + ml k 2 + m 2 J 1 ; 

k M* = k l_ ( AN ^f m 2 , K - <) 2 \ 

dk ~ 167r 2 a 2 \{k 2 + m 2 q ) 3 q {k 2 + m 2 n ) 2 {k 2 + m 2 a ) 2 J ' 1 ' 

As pointed out in the introduction, the inclusion of the wave function renormalization in 

the RG flow is crucial to obtain the correct behavior for the low /c-dynamics. We emphasize 

that the masses that enter on the right hand side of the flow equations (jMj) and (|35|) are set 

constant in the Mow and depend only on the parameter M 2 . This is indicated by using 

small letters for the masses. 

To compare this approach with chiral perturbation theory, we expand the pion decay con- 
stant, the chiral quark condensate and the pion mass in the current quark mass m c . We 
distinguish between the chiral symmetric and the explicitly broken system by keeping the 
current quark mass m c in all expressions. The parameter M 2 in the xPT-expansion is 
directly related to m c by the Gell-Mann-Oakes-Renner relation. In lowest order we identify 

M 2 = 2m c B, (36) 

where B is one of the low energy constants that appear in the chiral Lagrangian (JTJ. 
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For explicit symmetry breaking the effective potential £7($ 2 , a) depends on M 2 via the sigma 
field 



U($ 2 ,a) = U ($ 2 )-6a, 



A 



2^ 



f/ ($ 2 ) = - ($ 2 - $ 2 ) 



(37) 
(38) 



where we expanded the potential up to fourth order. We parametrize the mesonic fields and 
the symmetry breaking parameter as 



a 



71 



and S 



m 



uv 



guv 



(39) 



The chiral and physical minima of the effective potential are defined by 



and $ 



phys 




(40) 



which determines the pion decay constants in the chiral limit and the physical case as 



= V Zq, 0o and f n = V z $ a o- 



(41) 



We will denote the difference between the physical pion decay constant and the decay con- 
stant in the chiral limit by AL = f n — f n Q. The renormalized coupling constant is \r = -£r. 
Starting from the mass equations (J23|) . (|2*Bj) and (J27)l . we identify the M 2 -dependence of the 
quark and meson masses. For the pion mass we find 



m 



7r0 



1 d 2 U Q (<S> 2 ) 



m. 



1 d 2 U (<S> 2 ) 



tphys 



0. 



2X R f n0 AU + O(Af 2 ). 



(42) 
(43) 



As in chiral perturbation theory, the parameter M is the pion mass in lowest order of the 
chiral expansion 



ml = M 2 + 0(M 4 



(44) 



which enables us to identify 



M 2 = 2\ R U Af w . 



(45) 
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Analogously, we derive the expansions for the sigma mass and the constituent quark mass: 

1 d 2 U ($ 2 



m 



m 



o-o 

2 



Z^ da 2 
1 <9 2 f/ ($ 2 ) 



m 



qO 



m„ 



Z§ 

-— <& 



m'o + 3M 2 , 



^phys 



G 



G 



y phys 



Z q \JZ§ 



(f nQ + A/ ff ) = m g0 + A M 2 



(46) 
(47) 
(48) 
(49) 



where the coefficient of the M 2 -term of the constituent quark mass is given by 



A 



m q0 



fn g0 
mi 



(50) 



Once the values m a0 and m q o for the sigma mass and the constituent quark mass in the chiral 
limit are fixed, the M 2 -correction terms in m a and m q enter the flow equations. Although 
these terms do not contribute to the logarithmic corrections in m 2 and they influence 
the M 2 -corrections. 

We first discuss the results for the pion decay constant f n = o§\fZi&. To derive its flow 
we need the evolution equations for o"o and Zq>. With the parametrization of the effective 
potential [7($ 2 , a) (eq. (|3Tj) ). we can express the evolution of the expectation value through 
the evolution of the effective potential. The expectation value satisfies 

aw* 2 ) 



<T=o"o,7r 2 =0 



<9$ 2 



2a -6 = U^a 2 )2a -5 = 0, 



(51) 



since necessarily ft = 0. The prime denotes the derivative of Uo(§ 2 ) with respect to $ 2 . 
While the coefficient of the symmetry breaking term 5 does not evolve under the RG flow, 
all other coefficients in the potential Uo(§ 2 ) evolve, so that o~o = o~o(^) is a function of the 
RG scale. Details of the derivation can be found in appendix [X] Combining eq. (|5T|) with 
the flow equation for the derivative coupling Z$, we find 



k—f 2 



k 6 



+ 



+ 



1 



?7r 2 \(k 2 + m 2 ) 2 (k 2 + m 2 ) 2 (k 2 + m 2 )(k 2 + m 2 ) 

3 ml ( 1 1 
+ 



2m 2 \{k 2 + m 2 ) 2 (k 2 + m 2 ) 2 



--AN c N f m 2 



m 2 (k 2 + m 2 ) 2 



+ 



(k 2 + m 2 ) 3 



(52) 
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Integrating this equation, we find a logarithmic correction to f„ 



M 2 M 2 

U = U ( 1 - Y^glog^ + 0(M 2 ) ) . (53) 



The structure of the result is exactly the same as in x?T (cf. eq. (J7J)). Crucial for this result 
is the inclusion of the bosonic wavefunction renormalization, which enters at order 1/N C and 
is not recovered in the standard large Af c -limit. While the generic form of the logarithmic 
term arises entirely from IR pion fluctuations and is universal, the scale inside the logarithm 
as well as all other analytic corrections depend on the detailed dynamics at larger scales and 
cannot be trusted in this simplified analytic approximation. These contributions have to be 
studied in a detailed numerical analysis of the non-perturbative flow equations performed 
in the next section. 

Analogously to the pion decay constant, we can derive a flow equation for the square of 
the pion mass m\. As can be seen from equations and (|3*Tj) . the physical pion mass is 
derived from the potential U(§ 2 ), 

2 dU 



m 2 



*-* • (54) 



Differentiating the flow equation for the effective potential (j34j) with respect to $ 2 and 
combining it with the equation for Zq> (eq. (jSHJ) leads us to the flow equation for m\ 

d_ 2 dU'jal) 2 U\al) dZ* 2U"(a 2 ) da 2 

dk *~Z^ dk Z 2 dk + z„ K dk 

ml ( ( 1 3 m 2 \ ( 1 1 



167r 2 f 2 \\ 2 2 m 2 J \(k 2 + ml) 2 (k 2 + m 2 a ) 2 J (k 2 + ml)(k 2 + m 2 a 

AN r Nf 2 / 1 1 



+ (k 2 + m 2 g ) 2mq \k 2 + m 2 + m 2 J ) ' ^ 

We remind the reader that in this equation all values on the right hand side are considered 
constant in k and in lowest order in the chiral expansion, i.e. in particular m 2 = M 2 and 
f n = f n0 . Similarly to the case for integration gives us the logarithmic correction to m 2 

^4 + ^ ios £ +o(A 4 (56) 

As before for we also find for m 2 the correct logarithmic behavior compared with %PT. 
Again we do not calculate the contributions linear in M 2 , which we will compute numerically 
in the following section. To conclude, we emphasize that the RG approach to the NJL model 
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shows analytically the same low energy behavior as for the two observables f n and m 2 
considered here. It should be mentioned that the values of the zeroth order terms f n Q and 
M 2 depend on the specific values of the integration boundary k x sB- With the approximation 
of constant masses we cannot obtain reasonable zeroth order values with the same scale k x sB 
for /V and m 2 . For completeness, we have also derived the chiral expansion of the quark 
condensate (qq) in appendix El It also reproduces the logarithmic term known from %PT. 



V. NUMERICAL EVALUATION OF THE LOW ENERGY CONSTANTS 

In this section we present numerical results for the pion decay constant and the pion mass 
from the NJL model as functions of the current quark mass which parametrizes explicit chiral 
symmetry breaking. As demonstrated in the previous section, we reproduce the logarithmic 
dependence of / w and m 2 on the symmetry breaking parameter, since the fluctuations of 
the pions are correctly included. In chiral perturbation theory, the low energy coupling 
constants of the next-to-leading order chiral Lagrangian manifest themselves as scales A, 
in the logarithmic terms of the chiral expansion. These scales contain modifications due to 
additional terms proportional to ~ M 2 . 

For the numerical evaluation, we expand each coupling of the effective action (eq. (fl9|) ) in a 
Taylor series around the physical minimum $ p h ys = {&o{k), 0). As discussed in section UTTj we 
take for the wavefunction renormalization functions Z$ and Z q and the Yukawa coupling G 
only the first, ^-independent term into account. Since the effective potential U determines 
the vacuum expectation value a and therefore is the most important coupling, we expand 
U as follows: 

4 

U($ 2 ,a,k) = Y, a i( k )(® 2 -Mkfy -6<r (57) 

i=l 

where as usual $ 2 = a 2 + if 2 . The minimum condition for <Jo(k), 

= 0, (58) 



<T=CTo(fc) ,7? 2 =0 



gives as an additional constraint 



a ' w " (59) 

We solve the coupled set of four flow equations (|21| ) -()24| ) numerically for different values of 
the UV cutoff and for a wide range of current quark masses m c . To obtain its bosonized 
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form we perform a Hubbard- Stratonovich transformation of the initial NJL model. In this 
bosonized Lagrangian the kinetic term for the mesonic fields and higher meson interaction 
terms are zero, so that the meson potential is characterized only by muv- This gives the 
full set of initial conditions at the UV scale: 

l LrUV 

Z<s>,uv = 10 9 
G uv = 1 (60) 

Zq,UV = 1 

For a given cutoff Auv, we determine the parameter muv and the physical current quark 
mass m phys by requiring that the values f w = 92.4 ± 0.3 MeV and = 138.0 ±1.0 MeV are 
reproduced within these tolerances. This condition is the most natural choice and allows a 
calculation independent from xPT. As a consequence, the value of the physical current quark 
mass m p h ys varies with the UV cutoff Auv- To compare our results with chiral perturbation 
theory it is necessary to connect our symmetry breaking parameter m c with the parameter 
M 2 used in %PT. As already explained in section ITV] we set 

M 2 = 2m c B. 

In Figs. El and 01 we compare our results for f v and m 2 as functions of the current quark 
mass m c with those of x?T, obtained from j46j. We use a scaled current quark mass, since 
the absolute value of the quark mass does not have an immediate physical meaning and 
depends on a chosen scale. In ref. [46], the xPT results are plotted against the dimensionless 
variable m c /m s , where m c = m u = rrid is the mass of the u- or d-quark, and m s is the mass 
of the strange quark. The physical pion decay constant and pion mass are obtained for 
m c /m s = 1/26. 

In order to be able to compare our results with %PT, we introduce a scale m phys such 
that m c /m p h ys = 1 for the physical value f n = 92.4 MeV, and we rescale the x?T results 
accordingly. Over the considered range of UV cutoffs, the value of m p h ys varies from m p h ys = 
12.9 MeV for A uv = 1.0 GeV to m phys = 6.0 MeV for A uv = 1.5 GeV. By construction, 
the RG gives the physical value for f n at m c = m p h ys for all UV cutoffs Auv- As shown 
in Fig. El the curves L( mc ) become steeper with increasing cutoff Auv- For the cutoff 
value Auv — (1-00 ± 0.05) GeV, the RG result falls within the band given by %PT. In the 
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FIG. 2: The pion decay constant f n is shown as a function of the scaled current quark mass 
m c /mnhys for different values of the UV cutoff. The shaded region is given by xPT, obtained from 
ref. 46j]. The individual lines correspond to the full numerical RG results for different cutoffs from 
Auy = 1.5 GeV down to Auv = 1.0 GeV in steps of 0.1 GeV. The curve with the largest slope 
belongs to Auv = 1-5 GeV, the curve with the smallest slope corresponds to Auv = 1.0 GeV. 
All curves are fixed to agree at the physical point (7n c /m p hys = 1), which corresponds to f n = 
92.4 MeV. 

chiral limit m c — > 0, the pion decay constant from the RG and from xPT do not necessarily 
match. Depending on the value of the UV cutoff, our result for f^o = f n (m c — > 0) varies 
between 79 MeV for a cutoff of A uv = 1.5 GeV and 87 MeV for a cutoff of Auv = 1.0 GeV. 
In xPT, the chiral limit of the pion decay constant is = 86.2 ± 0.5 MeV jl7j ]. 
We note that the slope of as a function of m c j m p h ys increases with increasing value of 
Auv- in the NJL model, this is due to the increase of mjjy with the cutoff A^/y, which can 
be seen from Tab. HI As a result, for a given current quark mass m c , the linear symmetry 

Til? 1 

breaking term 5 = J JV m c also increases with the cutoff. This leads to the increase of the 
slope of fJ mc ) with A uv observed in Figure |21 

m phys 

The results for the pion mass are shown in Fig. El for the same set of UV cutoffs. In 
contrast to the pion decay constant, the slope of m^j. as a function of the symmetry breaking 
parameter decreases with increasing Auy. The RG result falls within the band of the 
result for a much larger range of cutoffs Auv, compared with the result of the pion decay 
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FIG. 3: The pion mass is shown as function of the scaled current quark mass m c /m p hy S for 
different UV cutoffs Auv- The shaded region is given by xPT, obtained from ref. The lines 

correspond to the RG results for different cutoff choices, from Auv = 1-5 GeV (smallest slope) 
down to 1.0 GeV (largest slope) in steps of 0.1 GeV. 

constant. Since the corrections to m 2 are already of order M 4 , and thus one order higher 
than the corrections to f n , m 2 increases almost linearly with m c /m p h ys . We find that for 
A uv « 1.0 GeV the results for /„. and m 2 are consistent with those of xPT, see Figs. |2]and 
|3] However, if one does not consider any additional physical observable, there is a priori no 
reason to favor any particular value of the UV cutoff for the NJL model, as long as the cutoff 
varies between 1 GeV < Auv < 1-5 GeV. A cutoff much smaller than 1 GeV is not justified, 
since the phenomenological A 4 -parameter is around 1 GeV. Consistent with this, ref. [47[ 
finds that for a cutoff slightly below Auv < 1 GeV the RG equations cannot reproduce the 
physical f n . A much larger cutoff would extend the effective four-fermion interaction into a 
region where dynamical gluon effects become important. 

While the leading-order behavior of f n Q and the slope of m 2 as a function of m c depend on 
the implementation of the symmetry breaking term in the model Lagrangian, the ratios 4 1 - 

JnO 

2 

and ^§ are more indicative of pion fluctuations (cf. Figs. HJ. They are determined by the 

pion mass in lowest order M and the low energy constants A3 and A4. 

The pion mass is normalized with M 2 = 2Bm c , where the value for 2B has been obtained 
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from the leading order fit to the pion mass. In the chiral limit, we have 

,2 



m c ^0 M 2 



1. 



(61) 



With this normalization, corrections to the linear m c -dependence can be compared with 
results of x?T j^j], see shaded regions in Figs.EJ In general, the cutoff sensitivity is enhanced 
in the expansion around m c = 0. At the maximum value of our calculation m c /m p h ys = 3.4 
it amounts to 20 %. 

In order to obtain the low energy constants, we fit the pion decay constant /^(f^c) and the 
pion mass mjirrir) (given in Figs. EJ) as functions of m c to the parametrization from ;^PT at 



one-loop order 



fir — fnO 

^ = 2B ( I 



2Bm c 2Bm c 2< 
log — — + 0{m Ci 



1 2Bm r 



log 



A 2 
2Bm c 

1V 3 



+ 0(m 2 c 



(62) 
(63) 



327T 2 fa 

We use current quark masses up to m c = 20 MeV, which corresponds to pion masses up to 
~ 250 MeV, depending on the exact value of the UV cutoff. 

In Table |U we show the pion decay constant in the chiral limit f n Q and the leading term of 
the pion mass M(m p h ys ) E 



/ 2fim p h ys as functions of the UV cutoff Auv, inuv an d m p hy S - 
As explained before, a fixed cutoff A uv leaves the two free parameters rxiuy and m phys , 
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FIG. 4: We plot the ratios f n / f^o and m^/M 2 as functions of the scaled current quark mass 
m c /m p h ys - The xPT results of ref. |4f| give the shaded regions, the RG flow equations yield 
the individual lines. For 7^/7^0, the RG results are shown from Ajjy = 1.5 GeV (top line) to 



A 



uv 



1.0 GeV (bottom line), and for m 2 /M 2 , from A uv = 1.0 GeV (top line) to A uv = 1.5 GeV 



(bottom line), in steps of 0.1 GeV. 



20 



Input [MeV] 


Output [MeV] 


Output 


ll UV 




' ' 'pnys 


f n 
JttO 


M(m u ) 




A„ 






i nnn 

1UUU 


933 9 


1 9 Q 


86 Q 


1 36 3 


43 3 


1 3D3 5 


9 99 


4 59 


1100 


265.3 


10.6 


85.1 


137.2 


147.2 


1297.6 


0.14 


4.49 


1200 


296.0 


9.0 


83.5 


138.5 


346.7 


1261.1 
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624.5 


1245.4 


2.99 


4.37 


1400 
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6.8 


80.7 
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931.0 


1247.8 


3.78 


4.37 


1500 


383.8 


6.0 


79.4 


141.0 


1217.0 


1261.3 


4.31 


4.38 



TABLE I: Values for the pion decay constant in the chiral limit f^Q, M(m p h ys ), the scale- 
independent low energy constants A3 and A4 and the scale-dependent quantities £3 and Z4 at scale 
M(m p h ys ) from the fits to the RG results, dependent on the UV values A[/y and m;yy. The values 
of rriphys are determined by the condition that our result should reproduce the physical values of 
/tt and m n at mph ys . 

which are then fixed from and The different values from Tab. HJfor different cutoffs 
lead to the following intervals for f^Q and M(m p h ys )- We find over the range of UV cutoffs 
considered 

79.4 MeV < U < 86.9 MeV, (64) 
136.3 MeV < M(m phys ) < 141.0 MeV. (65) 

Fits to our full numerical RG results with equations (JrT^ and (JHSJ) give the low energy 
constants A3 and A4, which are also listed in Tab. We extract from the table the ranges 

0.04 GeV < A 3 < 1.22 GeV, -2.29 <l 3 < 4.31, (66) 
1.25 GeV < A 4 < 1.30 GeV, 4.37 <l A < 4.52. (67) 

These ranges are obtained by varying the UV cutoff between the maximal and minimal 
values in the model, determined by the momentum region for which the NJL model can be 
considered as a suitable phenomenological description. Changing the cutoff, one finds from 
Table H] for example a large value for f n0 correlated with a small value of A 3 . The values of 
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the low energy constants used in x?T |46[ are 

0.08 GeV < Af T < 1.99 GeV, 0.5 < lf T < 5.3, 
1.12 GeV < Af T < 1.40 GeV, 4.2 < lf T < 4.6. 



(68) 
(69) 



The constant A 4 (Z 4 ) reflects the correction of the leading order result of /„. due to quantum 
fluctuations, and thus is less dependent on the particular implementation of chiral symmetry 
breaking than A 3 (Z 3 ), which determines the next-to- leading order term of the expansion of 
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FIG. 5: Comparison between the pion mass m n and the constituent quark mass m q , obtained from 
the RG flow equations, as functions of the scaled current quark mass m c / m phys- A large mass gap 
exists between m n and m q for small current quark masses ^ c / m phys- The plot is made exemplarily 
for a UV cutoff A uv = 1.0 GeV. 

Deviations between the results from xPT and the RG treatment of the bosonized NJL model 
are to some degree expected, in particular for large values of the pion mass. In the NJL 
model the low momentum regime contains free constituent quarks with masses of roughly 
~ 300 MeV (see Fig. HJ) for realistic values of the pion mass and the pion decay constant. 
In contrast to the assumption in xPT, the mass gap between the light pseudo-Goldstone 
boson m n and the more massive constituent quark m q becomes smaller for larger current 
quark masses, and the dominance of the pion field fluctuations is lost. This suggests that 
effects of quark loops or higher (qq) bound states may be enhanced for larger m c /m p h ys - 
Therefore we have investigated the range of m c /m phys up to approximately 13. For cutoffs 
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FIG. 6: The plots show the pion decay constant f n and the pion mass for large values of the 
scaled quark masses m c /m p hy S for the two cutoff choices Auv = 1.0 GeV(m p h ys = 12.9 MeV) and 
Auv = 1.5 GeV(m p h ys = 6.0 MeV). The shaded regions represent the results of xPT, while the 
bold curves give the result from RG flow equations. The curves are shown up to a quark mass 
m c ~ 170 MeV and m c ~ 80 MeV, respectively. Results from lattice calculations with Wilson 
fermions from |j are plotted for comparison and represented by the bold dots, while the grey dots 
give results from |5J] with staggered fermions. 

Auv = 1-0 GeV and Auv = 1-5 GeV, this amounts to current quark masses up to m c ~ 170 
MeV and m c ~ 80 MeV, respectively. The results for / T (ro c ) and m^im^) are shown in 
Figs. IH1 For A uv = 1.0 GeV, the expected variations for large quark masses do not leave 
the band predicted by For Auv = 1-5 GeV, the value of f w for large quark 

masses deviates more strongly from x?T. In the same figure, we show results of recent 
lattice results with Wilson fermions, taken from \m, which were obtained with a new pre- 

n 

conditioned Hybrid Monte-Carlo algorithm, and from |54J with staggered fermions. 
While the description of these lattice results does not require chiral logarithms, they do not 
rule out such terms, either. A comparison of large quark-mass lattice simulations with the 
curves of Fig. El allows an extrapolation of lattice data to small quark masses, as long as they 
lie in the region covered by our calculation and their slope is compatible with our results. 
A set of values for /„. and for large quark masses from our calculation can be found in 
table HTl in appendix iDl 
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VI. CONCLUSIONS 

We have studied the quark mass dependence of low energy observables within the N JL model, 
using renormalization group flow equations. We have shown that these equations analytically 
generate the non-analytic terms in the chiral expansion known from chiral perturbation the- 
ory. Our result confirms that the pion contributions to the low energy constants are treated 
correctly in the RG flow equations. For this it is essential to consider the renormalization of 
the kinetic term of the pion fields. While it is difficult in general to map non-perturbative 
RG results on a perturbative expansion, it is much easier to identify non-analytic parts. To 
derive these results, we treat the pion mass as a constant in the RG flow equation below 
the chiral symmetry breaking scale. Integrating the RG flow is then equivalent to perform- 
ing a perturbative one-loop calculation. The resulting chiral logarithms ~ M 2 log M 2 are 
accompanied by an analytic power series in M 2 , where M is the pion mass to lowest order. 
Going beyond the analytical results, we find that the numerical solution of the full set of the 
RG flow equations gives low energy constants which are compatible with the phenomeno- 
logical values from %PT. As input to these equations, in the bosonized NJL model one 
needs the triple of an UV cutoff Aj/y, a mass parameter muv and a current quark mass 
^phys- For a given value of the UV cutoff, the remaining input parameters are adjusted to 
reproduce the physical pion decay constant f n and the pion mass m 2 at the physical point. 
Therefore there is some freedom in choosing such a triple. Output of the RG approach are 
the values of the pion decay constant in the chiral limit f„o, the pion mass in lowest order 
M(m p h ys ) = A/2£>m,ph ys , and the two low energy constants A 3 and A 4 (see Table HJ). When 
the UV cutoff Ajjv is varied in a reasonable range 1.0 GeV < Kjjv < 1.5 GeV, the chiral 
limit of the pion decay constant varies as 79.4 MeV < f n o < 86.9 MeV and the pion mass 
in lowest order as 136.3 MeV < M(m p h ys ) < 141.0 MeV. The corresponding low energy 
constants are 0.04 GeV < A 3 < 1.22 GeV and 1.25 MeV < A 4 < 1.30 GeV. For the renor- 
malization scale-independent quantities we find —2.29 < h < 4.31 and 4.37 < h < 4.52 at 
the mass scale M(m p h ys ), see equation (JT3J). 

The NJL model lacks confinement, and therefore our model contains free quarks. As a con- 
sequence, the mass gap between the light pseudo-Goldstone bosons and the other hadronic 
states is reduced to a mass gap between the constituent quarks and pions. Nevertheless, 
our results show that the values of the low energy constants which we obtained from the 
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model are compatible with chiral perturbation theory. On the other hand it is well known 
that the linear a-model without quarks cannot reproduce chiral perturbation theory: Quark 
loops are essential to mimic effects of higher mass meson resonances like the p meson. 
The topic of low energy constants inphenomenological models has been addressed in many 
publication,, see e.g. Q Q ff U U fi Q Q. Ou r r esnl tS eo m p,e m en t an 
earlier study j3] that aimed at the low energy expansion within the linear cx-model and 
gave estimates for the low energy constants L 4 - L 8 in the 577(3) ® SU(3) parametrization, 
which correspond to li - / 3 , / 5 and l e in the S77(2) ® SU{2) case. These results showed a 
qualitative agreement with phenomenology, but relied on a mean field approximation and 
introduced an arbitrary matching scale. Since these assumptions are not necessary in our 
approach, it would be an interesting challenge to study also the low energy expansion of 
xPT at NLO. This can in principle be done without additional assumptions by an extension 
of our effective action to NNLO in the derivative expansion. 

The results presented here demonstrate that the NJL model gives a reasonable description 
of the infrared chiral dynamics and its quark mass dependence. By construction, it is based 
on the same representation of the chiral flavor symmetry as QCD. Therefore, it involves 
degrees of freedom that go beyond the low energy dynamics of the pions that are fixed 
by the symmetry. The requirement that a model must be consistent with the physical 
low energy constants provides constraints on the full set of possible dynamical theories at 
large momentum scales. We find that including quark effects, the results from the model 
are compatible with xPT, whereas a linear 0(iV)-Hiodel without quarks is not In this 
respect the NJL model passes an important test for applicability to chiral physics in the 
non-perturbative regime. This is also relevant for dense quark matter in the deconfined 
phase. Due to the difficulty of performing lattice gauge simulations in this regime, the NJL 
model is widely used as a tool for investigations of color-superconducting quark matter, see 
e.g. [^, and its phase structure, which is sensitive to the quark masses p\ . 
Finally, a quark mass expansion is essential to extrapolate lattice gauge simulations to the 
physical regime. The lowest-order expansion of chiral perturbation theory will be sufficient 
as long as the lattice masses are sufficiently small, but it ma y b ecome unreliable for large 
masses which are required for studies with chiral fermions j^, A description that 

explicitly includes quark dynamics beyond the infrared regime may give a better account 
of the large quark-mass behavior. We have given results for /„- and for large quark 
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masses in the curves in Fig. |U] and the table (Tab. |HJ) which can be used to extrapolate to 
the physical regime. The NJL model can help to bridge the gap between lattice simulations 
with heavy chiral fermions and the physical regime. 
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APPENDIX A: FLOW EQUATION FOR THE PION DECAY CONSTANT 

In the appendix, we have collected some technical details of the derivation of the various 
flow equations used in the approximate one-loop calculations. In this first section, we give 
the details for the derivation of the flow equation for the pion decay constant. We start 
by considering the minimum condition eq. (|51jl and take the derivative with respect to the 
renormalization scale. Note that 5 is independent of the RG scale: 

d 



k—u' 



4> 2 



2a Q + U'M) (k^a^j 2 a + U^a 2 ) 2 



k—U' 



2ao + 2^ [ k dk a ° 



0. 



(Al) 



We can identify the term multiplying the derivative of <7o as the mass of the sigma meson, 
evaluated at the expectation value ctq, 



Z*M'M) = 2U' Q {al) + 4o«Vo 



(A2) 



In this way, we can express the evolution of the expectation value in terms of the evolution 
of the derivative of the potential Uq(§ 2 ) as 



k—a 2 
dk 



4o$ 



Z$M 2 (a 2 ) dk 



(A3) 
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To get the desired differential equation for [7q($ 2 ), we differentiate the flow equation for the 
potential (J5%Jl and insert the derivatives of the masses from equations (|23j) - (j2*6*j) : 



k 6 f AN c N f dM 2 3 dMl 1 dMl 



32vr 2 [ (k 2 + M 2 ) 2 <9$ 2 (k 2 + M 2 ) 2 <9$ 2 (k 2 + M 2 ) 2 d<S> 2 

*• r mN f o> )f i,, (A4) 



32vr 2 i(A; 2 + M 2 ) 2 Z 2 2 $ 2 v CT ^ V( fc2 + ( fc2 + M «?) 2 , 

For the flow equation for the expectation value of the field we evaluate this equation at 
$ 2 = cr 2 . In order to obtain a flow equation for the pion decay constant, we now need to 
include the flow equation for the coupling Z§: 

k e { AN c Nj G 2 ji ( 1 2 \ 



16^ \~ (k 2 + M 2 ) 2 ^ \(k 2 + M 2 ) + M 2 J 
t M 2 a -Ml \ ( 1 1 
+ V 3 ^W| ) \{k 2 + Ml) 2 + (k 2 + Ml) 2 ' iA ' ,] 



(k 2 + Ml)(k 2 + Ml) 

APPENDIX B: THE CHIRAL CONDENSATE (qq) 

In principle it should be possible to derive the quark mass expansion of the chiral quark 
condensate (qq) analytically in a manner analogous to the pion decay constant f n and the 
pion mass m 2 in section HVI For this purpose one has to compute the chiral condensate from 
the partition function Z, 

Z[J, fj, rj\ = e"^ = J VqVqV® exp (-S[q, q,$) + j d 4 x( J$ + fjq + yfij , (Bl) 

by differentiating it with respect to the additional source term m c qq 

Using renormalization group methods, this would lead to a flow equation for (qq) similar 
to equations (J52|l and (J55j) for /„. and m 2 . Unfortunately, it turns out that we probably 
miss some relevant terms in equation (|B2|) . This may be related to a lack of rebosonization 
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of the evolving four-fermion term [59]. Therefore we restrict ourself to show that the (qq)- 
expansion follows from the other two expansions. We derive (qq) from our results for f n and 
m\ from section HVl by using the Gell-Mann-Oakes-Renner relation, = ^W-rn c , and then 

J TV 

replace m c by the lowest-order relation in the chiral expansion, m c = M 2 j^-. This gives 

/ ^M 2 M 2 \ 

= <55> »( 1 -32^ l0g *^ + ° (M2) )' <B3) 
As for the pion decay constant and the pion mass, we reproduce exactly the logarithmic term 
of xPT (cf. eq.(jHJ)). In the way employed here, this is a consistency check of our expansions 
for /,r and m 2 (eqs. (J53|) and with the Gell-Mann-Oakes-Renner relation. 



APPENDIX C: REGULARIZATION INDEPENDENCE OF THE CHIRAL 

LOGARITHMS 

In general, for finite values of the renormalization scale k, the results for the RG flow depend 
on the choice of regularization function. However, for the family of cutoff functions that we 
consider in this paper, 



{U [Tk2)lT ) 



f {a \rk 2 ) = [ Y^-{rk 2 yr \ exp(-rA; 2 ), (CI i 

we can show that the logarithmic contributions involving the pion mass are independent of 
the choice of function within this family. The general flow equation for the potential, which 
results from a cutoff function of this type for a > 2, is 

1 3 AN c N f 



kd k U k {&) = 

lb7T 2 a{a — 1) 



ik 2 + M2($ 2 ))«-! (k 2 + M 2 ^ 2 ))"- 1 {k 2 + M2(f 2 ))«-i 

(C2) 



As an example, consider the pion contributions to the flow equation for the expectation 
value of the field, <Jo(k): 

2Xa (k) 3 2 1 (k 2 ) a 
[kd k a oW \ pioD - mao{k)2) 1Qn2 * - a{k 2 + M 2 {ao{k) 2 )r ^) 

For simplicity we neglect the couplings beyond the four-point couplings as in equation 
and approximate ^Mg- = A. 
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The term relevant for the logarithmic contribution is 

1 (k 2 ) a 1 (k 2 + Ml- Ml) a 

a (k 2 + MIY ~ a (k 2 + MIY 

1 (k 2 + MlY - (?) M 2 (k 2 + MlY 1 + (|) (M 2 ) 2 (k 2 + M 2 Y~ 2 + - ■ ■ ■ 
" a (k 2 + M 2 Y 

- \ ~ M -WTW) + ^ {Mif WTMW " + " ' (C4) 

We note that the second term is independent of the parameter a. After setting Ml = M 2 
and on integration over the scale k, this term yields 

- l -M 2 \og{k 2 + M 2 ), (C5) 

and is thus responsible for the emergence of the chiral log. We conclude that in the current 
approach, the logarithmic contribution is independent of the regularization scheme. 
In general, the expansion of the observables in terms of the pion mass parameter M 2 is 
dependent on the choice of the cutoff function in our perturbative approximation. However, 
the leading logarithmic term in the pion mass parameter (the chiral log) is independent of 
the choice of a for the cutoff function. 



APPENDIX D: RESULTS FOR LARGE QUARK MASSES 

In the following, we list selected results for f n and from the RG calculations for large 
values of the current quark mass. We give results for several different UV cutoff values. 
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A uv = 1000 (m 


phys 12.9) 


A uv = 1200 (m 


phys 9.0) 


A uv = 1500 (m 


phys — 6.0) 


m c 


U 


m n 


U 


m n 


U 


m n 


10 


91.4 


120.9 


93.1 


144.8 


97.9 


175.9 


20 


94.7 


172.3 


99.3 


204.3 


108.0 


244.6 


40 


99.4 


247.1 


107.7 


289.1 


120.9 


341.1 


60 


102.9 


306.9 


113.5 


355.9 


129.7 


416.2 


80 


105.5 


359.6 


117.9 


414.2 


136.3 


481.0 


100 


107.5 


408.1 


121.3 


467.4 


141.4 


539.7 


120 


109.1 


454.2 


124.0 


517.4 


145.7 


594.4 



TABLE II: Values for the pion decay constant f n and the pion mass for large current quark 
masses m c from the RG calculation. The data are given for three different UV cutoff choices Aw 
All values are given in [MeV]. 
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